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We study gelation in suspensions of model colloidal particles with short-ranged attractive and 
long-ranged repulsive interactions by means of three-dimensional fluorescence confocal microscopy. 
At low packing fractions, particles form stable equilibrium clusters. Upon increasing the packing 
fraction the clusters grow in size and become increasingly anisotropic until finally associating into 
a fully connected network at gelation. We find a surprising order in the gel structure. Analysis of 
spatial and orientational correlations reveals that the gel is composed of dense chains of particles 
constructed from face-sharing tetrahedral clusters. Our findings imply that dynamical arrest occurs 
via cluster growth and association. 



Systems far from equilibrium often exhibit complex 
structures even while their equilibrium behavior may be 
quite mundane. A dramatic example of this phenomenon 
is the behavior of suspensions of attractive spherical par- 
ticles which at equilibrium phase separate but which, 
when concentrated rapidly, form a disordered arrested 
solid - a 'gel' - able to support a weak external stress 0. 
Although particle gels, formed from colloids or concen- 
trated protein solutions, play an important role in many 
areas of materials science as well as biology the molecu- 
lar mechanism of gelation remains far from understood 
0|- Recently, in an attempt to understand the origin 
of the arrested states which hinder protein crystalliza- 
tion, attention has focused 0, 0, S S S 13 on gelation 
in systems where short-range attractions (which favor 
particle aggregation) are complemented by weak long- 
ranged repulsions (which provide a stabilizing mechanism 
against gelation) . Experiments i^, i^, ^, and simu- 
lation [71 reveal rather unexpectedly a phase of stable, 
freely-diffusing clusters of particles. The appearance of 
an equilibrium cluster fluid has lead to suggestions 0, 13 
that gelation in these systems occurs via a cluster glass 
transition, where clusters (as opposed to particles) are 
trapped within repulsive cages generated by the long- 
range repulsion. 

In this letter, we explore experimentally the effect of 
weak long-range repulsive forces on dynamical arrest in 
a model system of uniform colloids with short-range at- 
tractions. We use three-dimensional fluorescence confo- 
cal microscopy to follow the process of gelation directly 
with single particle resolution. In sharp contrast with the 
structure of all gel networks studied to date, we find that 
the connected solid network exhibits an high degree of 
orientational ordering. Our findings suggest that kinetic 
arrest occurs via a one-dimensional cluster growth and 
percolation mechanism and not a cluster glass transition 
as has been proposed previously 0, 13 ■ 

Our particles are poly(methyl methacrylate) [PMMA] 
spheres stabilized against aggregation by a thin (~ lOnm) 
covalently-bound layer of poly(12-hydroxystearic acid). 
The spheres have a mean radius of a = 777 nm, a 
polydispersity of 3% about the mean, a mass density of 



Pc — 1.176 g cm^"^ and are labelled with the fluorescent 
dye DilCis 9]- The spheres are dispersed in a mixture of 
cycloheptyl bromide (CHB) and cis-decalin, which nearly 
matches both the mass density and the index of refraction 
of the colloid. In this mixture the PMMA spheres have 
a small but reproducible positive charge. Phase sensi- 
tive electrophoretic light scattering measurements reveal 
a net charge Q of -|-140e on each particle, equivalent to 
a surface charge density ~ 10^^ that of typical aqueous 
colloids. 

The potential between spheres, U{r), is the sum of 
long-ranged screened coulombic repulsions ULR{r) and a 
short-ranged attractive potential, UsR{r). The electro- 
static repulsion has the familiar Yukawa form, Ulr{t) = 
{A/r) exp[2Ka(l — r)]. Here the Debye length reflects 
the effective range of the coulombic repulsions screened 
by the ions in solution, r is the pair distance between 
particles scaled by the diameter 2a and A is the contact 
value of the potential, is considerably extended by 
comparison to aqueous systems as a result of the low 
solubility of ions. We estimate the ionic strength as 
around 5 x 10^^ mol dm^'^, corresponding to a Debye 
screening length = 1 ± 0.2/xm ^3- The potential 
at contact. A, was estimated from electrophoretic mea- 
surements as -fSOfcsT. To induce particle aggregation 
we generate a short-range attraction between particles 
by adding poly(styrene), a non-adsorbing polymer. We 
use a polymer with a mean molecular weight of 10^ Dal- 
tons and a mean radius of gyration rg « 92 nm. The 
attractive interactions induced occur at pair separations 
less than ^ — Vg/a or approximately 13% of the colloid 
diameter. Their strength is controlled by the polymer 
concentration and is equal to UsRlksT k, — 140p at con- 
tact, assuming the polymer behaves ideally; here 4>p is 
the fraction of the free volume of solution occupied by 
polymer coils. Again, we emphasize the markedly differ- 
ent range of the attractive and repulsive potentials; the 
attractions act on a length scale while the repulsions 
extend over a scale ko ^ ^a. 

To study the three-dimensional structure of particle 
gels we use fluorescence confocal microscopy. The sus- 
pension was contained in a cylindrical reservoir 10 mm in 



FIG. 1: Confocal microscope images of colloid-polymer mixtures at different volume fractions. From left to right: (j)c — 0.080, 
0.094, and 0.156. The attractive interactions are the same in all samples, Usr = — OksT {(j)p = 0.69). Images (a) and (b) 
contain clusters while (c) shows a network phase. The bars are 20 ^m long. Inset: Phases observed as a function of the volume 
fractions of colloid (^c) and polymer {(j)p). 



diameter and 650 fim deep created by sealing the edges 
of a 170 fim coverslip to a microscope slide. Fluores- 
cence was excited by a He-Ne laser (543 nm) and imaged 
through a 50 /xm pinhole. The digital image was opti- 
mized using the Nyquist criterion to select pixel size and 
frame spacing from the instrumental resolution. The mi- 
croscope's 63x oil immersion objective provided a field of 
view of 73 x 73 /itm^ with a magnification of 71 nm/pixel. 
Typically, 345 vertically spaced images were collected at 
a separation of 0.16 fiia. Three-dimensional particle co- 
ordinates were extracted with an accuracy of about 50 
nm using algorithms similar to those described in TTl . 

Confocal microscopy images of colloid-polymer mix- 
tures at fixed attraction Usb. = —9 k^T and three differ- 
ent particle concentrations are shown in Fig. ^ At the 
lowest volume fraction (Fig. ^) the sample contains a 
fluid of small, approximately spherical clusters. We con- 
firmed that the clusters were stable by repeatedly imag- 
ing the sample for over a week, observing no significant 
growth in the size of clusters. With increasing (pc the 
clusters grow in extent, becoming less spherical and more 
chain- like; the micrographs (Fig.^) collected 4 hours af- 
ter loading clearly show short strands of densely packed 
chains of particles with a relatively uniform thickness (~ 
2a). The orientation and center of mass of each strand 
fluctuates and the sample is ergodic. Finally, at high 
densities {(f>c > 0.1) the clusters form a connected solid 
network, which Fig. ^ reveals occurs with no change in 
local structure; the gel consists of linked densely-packed 
chains of particles, reminiscent of percolation. Although 
the long-range structure of the network is highly disor- 
dered, there is a surprising degree of uniformity at short 
scales. The chains which form the network are clearly of 
nearly equal thickness. Repeated scanning reveals that 
while individual chains display local fluctuations all long- 
range motion is suppressed and the sample is nonergodic. 
Using confocal visualization we followed the transition 
between the ergodic cluster fluid and nonergodic network 



as both the colloid density (f>c and the short-range attrac- 
tions Usr were varied. A well-defined boundary between 
the two phases is clearly evident (inset, Fig.0 which we 
identify as the gelation transition. 

The spatial correlations evident in the gel are intrigu- 
ing. The order is in sharp contrast to all previous reports 
of colloidal gels where the structure has been character- 
ized locally either as a scale-invariant fractal cluster (e.g. 
5]) or, less commonly, as an amorphous packing of hard 
spheres The confocal images confirm that the local 
order is not a consequence of gelation. The same local 
structure is clearly present in both the fully ergodic clus- 
ter fluid (near the gelation boundary) as well as the non- 
ergodic gel. Indeed the main change at gelation is the de- 
gree to which chains are interconnected rather than any 
change in their internal structure. Thus we hypothesize 
that a weak long-range repulsion coupled with a short- 
ranged attraction favors the formation of dense strands 
of particles, the length of which grow with increasing (j)c 
until, at the gelation transition, the strands become com- 
pletely interconnected and a gel forms. To test this hy- 
pothesis we determine the three-dimensional coordinates 
of our samples. We focus on the gels formed with attrac- 
tive interactions flxed at Usr = —13 k^T, referring to 
future publications for a more extensive analysis. 

We flrst provide numerical evidence of correlations 
that are independent of density and characteristic of the 
chains that form the network. We start by analyzing the 
translational order present. From the particle coordi- 
nates we compute the pair correlation function g(r) and 
hence the local density proflle (/)(r) = s'^g{s)ds, 
which is the average volume fraction within a sphere of 
radius r centered on each particle. Fig. HJa) shows (j){r) 
evaluated for different gel densities. The local volume 
fraction decays rapidly towards (pc on distances of 3-4 
particle diameters, indicating that translational order is 
short-ranged at all volume fractions. The situation is 
very different at shorter scales (r w 2a) where transla- 
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FIG. 2: Correlations in colloidal gels as a function of 
(a) The local volume fraction 4){r) and (b) the proportion o 
particles bonded to nt neighboring particles. Note that thi 
peak in ^(r) at r = 2a is independent of volume fraction an< 
the distribution of particles bonds does not change with Ac- 



tional order is enhanced; on these scales Fig. |2Ia) indi- 
cates that the local environment is dense and essentially 
independent of (f)c] the 'cage' around each particle has 
a volume fraction ^ 0.5 (from the peak of 4>{r)) and 
a size w 2a. To probe further the local organization 
we examine the orientational order. We define particles 
which are 'bonded' to each other by their separation r. 
If r < To, where tq is identified with the position of 
the first minimum in g(r), then the particles are con- 
sidered neighbors. This definition ensures that all par- 
ticles in the first coordination shell are counted as near 
neighbors. Fig. ^[h) shows the distribution in the num- 
ber of bonds per particle as the gel density is varied. 
The bond distribution is essentially independent of (/>c, 
confirming the picture of the gel as a network of invari- 
ant chains. The mean number of bonds per particle is 
rib = 5.6 ± 0.4, midway between the values of 2-3 ex- 
pected for a fractal structure and the 12 bonds found in 
dense liquids and glasses. More dramatic evidence for the 
three-dimensional order within the gel is shown in Fig.|3 
which shows a two-dimensional projection of particle cen- 
ters and bonds identified within a small slab of gel. The 
image reveals that the gel consists of relatively uniform 
chains constructed from dense clusters of particles. 

To check for local tetrahedral coordination we calcu- 
late the rotational invariants qi{i), wi{i) (for I = 4 and 
6), which are quantitative measures of the local bond 
topology around particle i. We define the bond-order 
parameters following 'l3'|. Structural and thermal fluc- 
tuations result in a distribution of 94, gg and wq. 
The resulting bond-order distributions provide a sensi- 
tive measure of the symmetry of the local environment 
|l3j . The wg-distribution, in particular, facilitates the 
detection of structures with tetrahedral symmetry since 
a large negative value signals a high degree of tetrahedral 
order. A tetrahedral cluster has wq = —0.106 while an 
icosahedron, an ordered arrangement of twenty tetrahe- 
dra, maximizes \'Wq\ at wq = —0.170 By contrast, 

the cubic symmetries found in /cc, hep and bcc crystals 
lead to near zero values of wq = —0.013, -0.012, and 0.013 
respectively. 

Figure 0] shows the bond-order histograms measured 




FIG. 3: A two-dimensional projection of the particle centers 
within a slab of gel (14.3 /im deep) at (j>c = 0.1. Particles are 
colored as a function of their depth within the sample and 
drawn 40% of their actual size for clarity. The bar is 20 /^m 
long. Inset: A spiral chain formed from tetrahedra of particles 
sharing faces. 



for the gel at 0c = 0.10 (the distributions obtained at 
other densities were similar but are not shown) . The pro- 
nounced maximum in the wg-distribution at ~ —0.13 
is particularly striking and demonstrates a high degree of 
orientational order. The peak is at too low a value how- 
ever to be accounted for by tetrahedral units alone and 
the small number of bonds per particle is incompatible 
with an icosahedral symmetry. Although neither struc- 
ture is consistent with the data, the position of the wq 
maximum, intermediate between the values expected for 
tetrahedral and icosahedral clusters, suggests that the gel 
probably contains ordered clusters of tetrahedra. 

Polytetrahedral order has been invoked to understand 
the structure of quasicrystals, atomic liquids and glasses 
[15I UtI however the concept has not been applied 
to gels, as far as we know. A viable structure must ac- 
count for the chains evident in the confocal images to- 
gether with the strong propensity for six-fold coordina- 
tion found experimentally. These constraints eliminate 
many possibilities. So for example, while linear chains 
of edge- or corner-sharing tetrahedra are common motifs 
in structural chemistry the number of bonds per particle 
varies along the chain and averages out at less than 6. 
More dense tetrahedral chains can be built by sharing 
faces of neighboring tetrahedral units. There is a sin- 
gle dimer of two face-sharing tetrahedra with the sym- 
metry of a trigonal bipyramid, one trimer (a bicapped 
tetrahedron), two distinct tetramers and five pentamers. 
The number of such arrangements grows rapidly as more 
tetrahedral units are added so that an infinite sequence of 
face-sharing tetrahedra may exist in any one of an infinite 
number of possible configurations, each with a different 
chain contour. Each configuration however has an iden- 
tical shell of six nearest neighbors since the neighbors of 
each particle belong to a trimer of adjacent tetrahedra, 
which may be realized in only one way. Here we focus 
on the most striking configuration, the graceful Bernal 
[3 (or Boerdijk-Coxeter [llll^) spiral, reproduced in 
Fig. 13 We have calculated the distributions of bond- 
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FIG. 4: Orientational order in a colloidal gel with <j)c = 0.10 and Usr = — 13k_gT. Measured wa, wa, qe and g4 bond-order 
parameter distributions (points) together with values calculated from a model of face-sharing tetrahedra (black curves). 



order parameters from computer-generated coordinates 
for a spiral strand consisting of 1000 tetrahedra. Ther- 
mal and structural fluctuations were simulated by adding 
random deviations, chosen from a Gaussian distribution 
of width Aa, to the ideal coordinates. Fig. 01 depicts the 
bond-order distributions calculated for the Bernal spi- 
ral with A = 0.21 where it is clear the measured wq- 
and W4-histograms are reproduced rather well by this 
model, strongly supporting our picture of the gel as a 
network constructed from chains of face-sharing tetra- 
hedra. The measured qe- and g4-distributions are less 
well reproduced but this is probably linked to the fact 
that our model only has six-fold coordinated species while 
the measurements contain particles with both more and 
fewer than 6 nearest neighbors. 



In summary, we have shown that the mechanism of 
dynamical arrest is profoundly altered by the presence 
of long-range repulsive interactions. Three-dimensional 
confocal microscopy reveals the formation of stable chain- 
like clusters, constructed from face-sharing tetrahedral 
units, which associate into a fully connected network 
at gelation. The possibility that a delicate balance of 
attraction and long-range repulsion may stabilize equi- 
librium clusters has been discussed in recent theoretical 
work ji^, 0| . However it is not clear if the formation of 
tetrahedra reflects the presence of charge or is simply a 
consequence of a centro-symmetric attractive potential. 
An explanation for the unexpected structures seen in our 
experiments remains an open challenge to theorists. Fi- 
nally it is tempting to draw analogies between the one- 
dimensional growth seen here and the tendency to chain 
formation observed in several protein solutions |U . 
These similarities suggest that weakly-charged colloidal 
systems could offer valuable new insights into the micro- 
scopic mechanism of protein nucleation and gelation. 



We thank Macolm Faers, Adele Donovan and Laura 
Starrs for extensive discussions and help during the ex- 
periments. This work was supported by EPSRC, Bayer 
Cropsciences and MCRTN-CT-2003-504712. 



[1 

[2; 

[4: 

[6: 
[7: 

[8 

[9: 

[10 



[11] 

[12 
[13 

[14 

[15 

[16 
[17 

[18 
[19 

[2o; 

[21 
[22; 
[23; 



K. A. Dawson, Curr. Opin. Colloid Interf. Sci. 7, 218 

(2002) . 

See, e.g. Non-Equilibium Behavior of Colloidal Disper- 
sions, Edinburgh, UK, 2002, [Faraday Discuss. 123 

(2003) . 

A. Stradner, H. Sedgwick, F. Cardinaux, W. C. K. Poon, 
S. Egelhaaf, and P. Schurtenberger, Nature 432, 492 

(2004) . 

H. Sedgwick, S. Egelhaaf, and W. C. K. Poon, J. Phys. 
Condens. Matter 16, S4913 (2004). 

A. D. Dinsmore and D. A. Weitz, J. Phys. Condens. Mat- 
ter 14, 7581 (2002). 

P. N. Segre, V. Prasad, A. B. Schofield, and D. A. Weitz, 
Phys. Rev. Lett. 86, 6042 (2001). 

F. Sciortino, S. Mossa, E. Zaccarelli, and P. Tartaglia, 
Phys. Rev. Lett. 93, 055701 (2004). 

K. Kroy, M. E. Cates, and W. C. K. Poon, Phys. Rev. 
Lett. 92, 148302 (2003). 

A. I. Campbell and P. Bartlett, J. Colloid Interface Sci. 
256, 325 (2002). 

The ionic strength was determined from the measured 
sample conductivity (170 pS cm~^) assuming ion diffu- 
sivities of 6 and 2 xlO^^" w? s^^, typical of Br^ and 
alkyl cation diffusion. The Debye length was calculated 
at 4>c = 0.1, allowing for the screening from the solvent 
ions and the Q counter-ions associated with each particle. 
J. C. Crocker and D. G. Grier, J. Colloid Interface Sci. 
179, 298 (1996). 

P. Varadan and M. J. Solomon, Langmuir 19, 509 (2003). 
P. R. ten Wolde, M. J. Ruiz-Montero, and D. Frenkel, J. 
Chem. Phys. 104, 9932 (1996). 

P. J. Steinhardt, D. R. Nelson, and M. Ronchetti, Phys. 
Rev. B 28, 784 (1983). 

J. D. Bernal, Proc. Royal Soc. (London) A 280, 299 
(1964). 

D. R. Nelson, Solid State Phys. 42, 1 (1989). 

J. P. K. Doye and D. J. Wales, Phys. Rev. Lett. 86, 5719 

(2001). 

A. H. Boerdijk, Philips Res. Bull. 7, 303 (1952). 
H. S. M. Coxeter, Can. Math. Bull. 28, 385 (1985). 
S. Mossa, F. Sciortino, P. Tartaglia, and E. Zaccarelli, 
Langmuir 20, 10756 (2004). 

J. Groenewold and W. K. Kegel, J. Phys. Chem. B 105, 
11702 (2001). 

M. Weijers, R. W. Visscher, and T. Nicolai, Macromol. 
35, 4753 (2002). 

C. Le Bon, T. Nicolai, and D. Durand, Int. J. Food Sci. 



and Technology 34, 451 (1999). 



